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ABSRACT 

Ali present stiffened cylindrical shell design formulas 
Memmedc Cesc Ol exvcinas hyGrostatic pressure were surveyed. 
ieeeweiso incluaed present desigcn practices for allowances 
gue to pressure hull imperfections, and actual test dat 
when available. Formulations for all.three basic null ° 
failure modes were then selected, Piroueeor Sccurady, and 
Secondly for conpatiblility with both elastic-perfectly 
Pieepic and stralin-hardenins met tals, whenever possible. 

tae 2ormulas were tnen inserted in & Logical flow 
pattern to design PlgetC- Joy, i ceuly plastic scantlings for 
failure ot the most efficient collanse mode. The process 
was programmed in FORTRAN IV and designed to iterate, varying 
several scantling verameters Be ti eatin. The 


"optinun" design, based on &@ simole hull weight/bduoyancy 
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PeGio Macmoelected at the comoletion of the run. 

Peeetoam inputs are: collapse depth, hull diameter, 
hull length, internal bulkhead svacine (svecified or 
unsvecified), framing (internal or external), and metal 
properties. Outputs are: shell thickness, typical frame 
spacing, typical frame size, heavy frame (bulkhead) spacings, 
and heavy frame size for each design, and an optimum design 
designation. Simple directions are given for conversion 


of the program to one compatible with strain-hardening métals, 
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INTRODUCTION 


one 


The amount of literature concerning tne collapse of ring- 
Stifiened cylindrical shells under hydrostatic pressure 
accumulated in the last fifty years is voluminous. This is 
understandable; the subject is very involved. To this day, 
Pex ect SOlUtiOon for all aspects of the problem does not 
exist. Good solutions for the different failure modes do 
earect, tLhougn, modified in varying amounts by empirical data. 
No attempt will be made here to list or summarize this 
Knowledge. Many have elready done this, and the finest 
review to date in this author's knowledge has been done by 
Mie G. Pulos?!® for the Navy's former David Taylor Mhodel Basin. 

It appeared that one should be able to integrate Hen 
albeit incomplete, yet extensive knowledge with the use of 
present generation computer science. Hand calculations for 
only one seometrical combination of shell thicsness, frame 
size, frame spacing, hull diameter, hull length, etc. are 
notoriously laborious, even for only one mode Ory ©24. re. 
submarine design processes using the hand technique would 
acnieve Beciate Sere res, DUS With Little or no idea if 
anything better existed. : Optimization, with the excevtion of 
a few conbinations tried at a Picea. Cosc Jn time, was ouv of 
mee scvuesv1On: wiile Similar Submarines covld be designed on 
past Knowledge, dirferent hull geometries or deeper operating 


depths meant a great deal of time and work, It was at the 





suggestion of Professor Evans that the development of a 
computerized Gesien ODViImizZetion was undertaken. The basic 
@eeien coueations were there. The computing tools were 
readily available. All that needed to be done was an 
[/ivesracvion of the two into a practical, useable, and most of 
all, reliable (in terms of latest empirical data, if 
necessary) program. 

EMoruly afvuer une start of the project, it was 
discovered that a smiliar progrém nad just been completed’? 
for the Naval Ships Peeearch and Development Center (NSHDC). 
ie. 1S hoped that by using some different approaches and 
techniaues, this program could be a valuable tool to use in 
Pemejicci0on with reference 22. 

~ The description of the program development will be done 
by SUudprograns, each building on the other, and ending with 
the optimization scheme of the main program. While the 
program develoved can be used only with elastic-perfectly 
Bolestic, isotrovic materials (e.e., HY-80 steel), it can 
easily be modified so as to be apvlicable also to strain- 
merdening metals, such as HY-150 steel. Furthes discussion 
on this will foliow later. Included also in the thesis will 
pe Various data obtained usine mie Oroeram or porvions 


Pmereo: ln parametric-type studies. 





It is recommended that wnile reading throush the various 
@eepecrs On SUbprograms, reference be made to figure 23 
Memeover 10), which is a simplified main flow diagram for 


the entire optimization. 
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CHAPTER ONE 
BACKGROUND INFORMATION AND TERMINOLOGY 
The entire hull design program is based on the three 
fundamental failure modes for ieee erence eyvlindrical 
meeieec. 2nerefore, in order tnat terminology remain 
consistently clear throusgnout the discussion, a brief 
description and catagcrization of these modes follows. Refer 
moriipure one for pictorial representations, 
fe AXLSYUMEDRIC FAILURES BETWEED STIFFENERS 
All axisymmetric failures, whether elastic, plastic, 
or some mixture thereof, are characterized by one or more 
Meecoruian--like pleavs, or circumferential ripples between 
maememrranes, For true axisymmetric failures, the stiffeners 
eosin undeformed. 
ee om eee ELD) PRILURE. This type of failure 
CCecUrocmomemyi1 bn elastic-perfectly plastic 
(plateau-type stress-strain curve) materials. It 
is regarded, then, as almost totally a yleld-type 
failure, although it is initiated parvo ly by 


tabilityvy phenomena, 


{unto 


ag) 


mn 


Cec ome Celie STC FATLURE. Inelastic 
failure is Also a failure ebove the purely 
G@lastic range, but in strain~nardening materials: 


Thus, the failure is in the range where Young's 


Mecano yes. and can intuitively be 





considered a kind cof combination elastic-plastic 
failure. Often this failure can be at a lower 
pressure than a pure buckling (elastic) failure, 
due to certain combinations of low modulus and 
hull geometry. The pure yield failure (A.1.) 
cannot occur in stra@in-hardening materials, 

pee ow ethic ELASTIC FATLURE. This is 
axisymmetric failure in the Lines ely elastic 
PonecwmmmeenicOrevically, this could occur, .given 
Peoveuerconerry, Ln eitner of the two above types 
of materials. Generally speaking, the required 
BeComeumeapmamone Of & thin shell relative to the 
NUL ete ter and deptn. in reality, this 
failure is, at this writing, a mathematical 
phenomenon only. Other hull failure bees Coc, 
rirst, Someacs ees has never been acnieved in : 
actual testing Tite is valuable, howevezx, in 
Getermining effects of geometrical defects on 
collapse pressure. 


B. ASYMMETRIC (LCBAR) BUCKLING FAILURES BEUVERN 
SIJFFENERS: 


HL lobar buckling failures are characterized by 
lobes of buckling distributed partly or completely around the 
shell circumference. The frames remain intact. 


eee emo eC HODES: This is primarily a 





DVvGr ees sal lure. The collavse pressure is 
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dictated by the modulus (elastic or inelastic) 
and hull geometry. 


C. GENERAL INSTABILITY (SHELL AND FPANE) FAILURE 


site ty doa tod st EN dls sai PEE EE iE ay a Pipe pt cid 


ts 


General instability failures are characterized by 
failure of both oe and shell simultaneously, sometimes 
extending the entire lengtn of the cylinder. 
ee eemo es INE STIG HODES. ‘hese are also 

bMemempemicl lures, except that lobes extend 
Lone budinally as well as Cin elimeerentialily, 

It is normally assumed tnat only a one~half wave 
extends between heavy frames (oedmcads ).. since 
these heavy frames are designed heavy enough not 
to deform at the general instability collapse 
DreSSurecmmeOLS combinavion failure of frames 

and shell is not as well understood ocr formulated 


mathematically as the other two modes, espneciaily 


v 


- 


In tUnewimelaesuic resion, Thus, larger safety 
faACtors @fe employea wnen checking hull designs 
in this mode. 

It should be emnhasized at the outset that the program 
to be described designs the stiffened cylindrical snells to 
Boerein the axisymmetric yield failure mode, For the 
Mntvested inelastic portion of the program, the design failvre 


mode would be axisymmetric inelastic failure AW ches 


PMilosoony 1S used vocey in submarine design end is 
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edavocated by most naval arcnietects. It has two basic 


PEaASONS: 


imme eminre by yieldine of the shell utilizes full 


2) 


material (i.e., yield) strength. Stresses in 
buckling type failures are usvally below yield 
stress, with failure depending mainly on the valve 
of Young's modulus and geometry of the failing 
structure. More efficient structures should thus 
Peewee rom Gesigning to &@ yielding failure. 
Imperfections in construction (e.g., hull out-of- 
roundness) effect buckling failures far more 
seriously than yield tines. = Haier rut Nes 
designed for yield failure would have less 
stringenv requirements for building, and less 
CnancewO micitiure, piven that imverfections 
might exist. 

TiewowitmerZ2cvlOn Criterion used is a simple 
hull weight Wem oul! sdgeplacement ratio, wnich 
appears to be the best general measure of design 


Cu le@e lemme ois S type of sS®ructure. 


e 
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Pee Oat Oli PENED CYLINDRICAL SHELLS 
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AXISYMMETRIC YIELDING (BETWEEN FRAMES) 


COLLAPSE FOEMNULAS: LUNCHIK(DTINB# 1291,1393) 





ASYMM®TRIC(LOBAR ) BUCKLING (BETWEI EN FRAMES) 
COLLAPSE FORMULA: REYNOLDS(DTNB#1 392 ) 





Chie eee eter ty (CONCURRENT FAILUFE 
OF SHELDE AND 2 


COLLAPSE FORMULA: KRENZKE-hIERNAN(DIME# 1677) 
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Cierink TWO 
peo oot 

Numerous subprograms througnout the optimization require 
tne values of various stresses in the hull. These stresses 
meeruae Circumfercntial and longitudinal stresses, at the 
Premes and av midbay, at the inner and outer shell surfaces. 
Several solutions for stresses due to external hydrostatic 
pressure on ring-stiffened cylinders have ayer eee ii ee 
past. The most famous was due to the Germans, Von Sanden 
Pema vnvner, in 1920. Portions of their analysis are still 
mamise vLoday. in 1930, the Italian Viterbo modified their 
analysis to include the so-called stiffener-expansion effect 
Meresult of axial stresses in the shell). Neither of 
these carly analyses, however, included the "bean column" 
meeect. This effect, introduced by Salerno and Erie. in 
1951, is caused by the interaction of longitudinal bending 
and lonsitudina) compression in the hvll caused by tne axial 
Mometo. Of the hydrostatic pressure acting on it. ‘The 
Salerno-Pulos stresses are an exact solution, and tne bean- 
column effect BCcCOUunvS for any non-linearities between. 
pressure and strain in the cylinders. In all cases except 


one (see chapter four), the program uses the more accurate 


Balerno-Pulos (hereafter S-P) stresses. 
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The beam-column effect is represented in the S-P 


analysis by the parameterV (nereafter: GAMA), where: 





Zz 
PeteP. O ~ (Se) 
Sp cevo-Y? \4 


16,18 oe t? . ° a) a) 
e! Ssmpne Terres cal )0ag cor 


Pee Oct ine 
Peemayimetvric elastic buckling of an unstiffened cylindrical 
cert under the action of uniform axial pressure." GAHMA=0 
corresponds to a zero beam-column effect (i.e., the stress 
solution of Von Sanden and Gunther, hereafter V-G). As 
GAMA grows larger, the beam-column effect, and thus the non- 
linearity between pressure and hull stresses, increases, 

When GANAS 1, then theoretically the between-frame failure 
mode snifts to axisymmetric elastic buckling (sce chapter 

One for definition). As explained in Chapter i, this appears 
to be a2 mathematical phenomenon only, for it has never been 
achieved in actuel testing, Other modes of failure (e.g., 
lobar buckling) always appear to occur first, or else the 


a4 


axisymmetric failure is always accompained by some yielding: 
Wem Stresses ere calculated by placing various 
combinations of S-P "F functions" (see reference 18 for 
expressions and curves) into Dae SP stress equations. 
Formerly, tnis was a very laborious process, far more time- 
consuming than obtaining stresses tnrough the simpler V-G 


equations. Tne curves developed by H.4. Krenzke and R. D. 


v 


Bhort (included in reference 18) snortened the labor 





N 
aS 


@ensiaerably. The computerized solution makes it almost 
IMandatory to use tne sunerior S-P stress solution, 

The "F functions" are obtained through the use of 
two basic subprorrams (see figures 2, 3 and 4). Subroutine 
FRAME essentially takes snell thickness, unsuopvorted shell 
length, and frame moment of inertia and manipulates them to 
Obtain the S-P variables THETA, ALFA, and BETA (see Appendix 
A for method of obtgining frame dimensions). These variables 
meemcnen transferred to the stress program PULOS via a 
COMMON statenen Pe oom secparacve 1rom PULOS because in 
one program, (HVYFRM) PULOS is used in two different places 
meme tne same scantlingss, All variables in FRAME are 
computed in accordance with the S-P stress analysis” ~, with 
the exception of the expression for effective een ee 
@eeerence 18 lists this as (in prograre terms): 

AR=AP2(EM/EF), 
AN=AFPs(EM/HF)se#2, depending on whether 
Memrramine is internal or external, resvectively. This is 
mo | strictly correct, and has since veen refined by Shore | 
wero used the sinilar equation: 
AG=APS(EM/PRF )eeQ, where Q=1+ 2GlhU 

Boies equation is good for eitner internaily or externally 
freaned cylinders with "reasonable" (i.e., suitable for this 
program's vurposes) ee Geoun/shiell radius ratios. 


Mmemsvoroutine PULOS is & straightforward acavation 


@emeune o-Y stress analysis. It produces the four "F 
Meret ions" and the variable "a" (see reference 16), the 
Various combinations of which are used in several parent 
prosrams to compute the hull stresses. 

irom the expression for GAMA, it can be seen tnat an 
Memeo pressure is required for PULOS to compute its outputs. 
imeotvner words, FEAMNE and PULOS can be used to calculate 
Mie Null stresses, given scantlings, material properties, 
mmo nydrostatic pressur more Of len, however, the progran 
Heeaeycvemoting to find a critical failure pressure. This 
resvits in two unknowns, the pressure and the stresses. 
GAKA thus becomes the result of a transcendental process, 
in whicn PULOS calculates stresses, and the pressure input 
is from a parent program dealing with a particular null 
failure mode, 

Geeay Cifficulvies Were experienced in this iterative 
process when GAMASi.0; this meant that ETAL (see figures 
3 or 4) became imaginary (i.e., the failure mode had shifted 
bo axisymmctric elastic). To combat this, several methods 
were employed, each one being different ror; dir lerent 
Meee. prosrans, This is the reason for PULCS and PULOS li. 
Maeemeotvhnods of validly circumventing this pitfall wiil be 


Peolained separately under each individual failure mode 


POP PrOsTran. sa 
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CHAPTER THREE 
COLLAPSE PRESSURE REDUCTION 

Meatasvacelly speaking, no cylinder that is manufactured 
today can be considered "perfect". There will always be 
some reduction in the collapse pressure Gue to manufacturing 
imperfections, such as shell or frame out-of-roundness, and 
to residual stresses from welding. If, however, a pressure 
hull is machined rather than rolled and welded, a structure 
Meee is “perfect" for all mie ties | purposes may be attaincd. 
Of course, for large pressure hulls, the expense (or even 
impossibility) incurred due to size prohibits construction 
by machining only. Thus, some allowances must be made in 
scantling computetion. | 

The best overall method (i.e., including perfectly 
Bleastic and strain hardening materials) yet devised is 
presented as a graph in reference 10 (see figure 5). Here, 
the lower curve represents an envelope of numerous model 
test resuits conducted over the years at the ilodel Basin. 
it may be observed that "the factor wnich lis a1) important 
in determining imverfection sensitivity is the margin of 
stability Do/5;, the abscissa on the plot. The lower the 
margin of stability, the greater the sensitivity to 
imperfections.“ The reduction factor (FEDFAC, OF Pexe/ Pexp 


for machined models }) located on the ordinate is applied to 
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collapse pressures of all modes in the optinization 
program (see 

the envelope curve is approximated by the equation: 

TEDPAC=0,6466740.1136725T0-0. 90096742 TO##2, where 

RTO=pe /D4 

This equation and the general form of REDPR were adopted 
from a similar program in reference 22. 

meommentioned above, this safety reduction factor is 
meemicable not only to this optimization, but also to one 
for strain hardening metals as well. For this reason (i.e 
miembect it is easily convertible to the strain hardening 
case) it was chosen above other existing wiliceiPo ome 
@nalyses for the strictly elastic buckling cases. No 
mathematical analysis presently exists for the inelastic 
case, due to its complexity. 21 Purener Giscuss10n on this 


memcontained in Chapter 13. 
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CEAPTAR FOUR 


ASYMNETRIC (LORAR) BUCKLING 


The most widely used evans ere for this failure mode 
eee evuie years has been the so-called "DINS Instability 
Formula", developed by Windenburg and Trilling. The formula 
was good only for an elastically-perfectly plastic material. 
More recently, Reynolds”? has developed a more generalized 
formulation wnich may be used for either elastic-perfectly 
plastic or strain~hardening materials. In reference 20, 
Reynolds recommended using the V-G stresses (longitudinal 
and circumferential at mid-bay). He stated that the accuracy 
of the analysis was not seriously impaired by not using the 
more Peper sone, yet more accurate S-P stresses. This is 
Peeaonice lly borne out by the exemple given in fissure 7. 

In any case, the theory correlates very closely with 
experiment, to within four Meee. 

The S-P stress programs (FRAHE and PULOS) can easily be 
made common with any other subprogran, Origsinellv it Was 
decided that there tiovuld be redundancy involved if the less 
accurate V-G stresses were used for the Reynolds lobar 
buckling subvrogran (HNLDS). Difficulties, however, were 
immediately encountered wnen tne S-P stresses were utilized. 


‘ 


Pres is used in the main program chiefly os a checking 


0 IO IE 


MeaCuLaie oince the hull scantlings, as mentioned in 
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inapter one, are designed for yield failure, ENLDS is used 
eory to insure thet the scantling set under scrutiny does 
meeeta2l by lobar buckling. “hus, since scantling sets 
being examined will usually fail firey ey sax Ley carie Uiclc 
yielding, failure by lobar bucklins sometimes occurs at much 
greater pressures. In many cases, failure by axisymmetric 


elastic buckling will occur at a pressure betuecn failure by 


{ 


Meld and failure by lobar buckiing. £s noted in chapter 
mo, Lhis phenomenon causes GAIL’ (beam-column effect) to 

increase in value over 1.0, resultins in imaginary values 
occuring within the PULCS SUDDLOSTAMS « PieGorts COL Torce!| 


convergence of HNLDS by inserting, for instance, values of 
GAMA=0.95, or of taking only absolute values of the eeGher tl 

Pieo-GAlit, were only mildly successful. The value of final 
convergence in @ny case was not accurate, and certainly was 
not that of lobar collapse pressure. 

For these reasons, then, V-G stresses were used in 
PNLDS (see fisure 8), as originally recommended in reference 
20. One distinct advantase of the V-G stresses is that 
mievercoulre no iteration for convergence. There 1s no 
separate stress DIOKTAN needed, and the V-G stresses ®re 
directly (and quickly) computed within HNLDS. Since (see 
chapter ten) RALDS is used itself in an iterative process 


Within the main prosram, this results in suovstantial savings 





O4/TB/7C ole, 


AUICFLCh CHART SET - FURKCTION RNLES pace ~ oF 


CHART TITLE - FUNCTICN RFNLOSET,SL) 


CCMPUTES CRITICAL 
ASYFMETRIC [LCBAR) 
BUCKLING PES SURE 


Re = &'/2.0 


4S 


SF & 2.7°5C6-Oer? 


enen ©? ep ap oe oe ae ae 


J 

] 

I 

J 

Y= 1 
Ro@BecG@°rqomety f 
J 

j 

! 

! 

! 


Foe & 0 8 


® 
' 
4 
' 
t 
‘ 
’ 
4 
’ 
‘ 
4 
' 
é 
4 
4 
‘ 
4 
4 
4 
' 
* 


2 
' 
' 
| 
s 
' 
t 
d 
4 
4 
U 
I 
] 
‘ 
‘ 
é 
i 
‘ 
4 
' 
td 


CaAu%92)) 
SSC eC E*SL/USCRT 
(ROT) ) 


STNHOx,¥) 
Yi/2. 


o 


COSHE X,Y} 
Y1/2.0 


we we we wm ew eee eww ee 


n 
~ 
al 

e 


1 
- 
- 
' 
me es we ws = ee = = ce 


TP2 = 
EXPITETA/S2 .0} 


IN2 = 1.0/TP2 


BEETAL = ! 
TETA/S2.080(S TAH $ 
{1P2,7TN2) ¢ | 
SUN(TETA/S2.09) | 
SECCSH(TP2,TA2) | 
COSETETA/2.C))) | 
J 

| 

| 

] 

* 


IP = TP2*#2 


Pim @ BF wee = oe = * 


ee wm ee ee oe ee 


| BEETA = ] 
§) TETA*LESINHE TP, J 
| IN) ¢ | 
| SINCTETADD I 
§ “SECOSHETP,TNI ~ | 
I CCSETETAI)?) J 
s 


V¥o=) Ol O/ (lee 
{1.0 - | 
GNL/2.0) | 
*CAF/ (SL 4T)) | 
#(BHETAL ~ j 
OoS®RHETAI/E1.0 + | 
C.5tSRETARCAF © § | 
Bet /¢Steti)) | 


me we we we we we a a ee eee 


mee ee gs es = 


I PCRE = 1 
§ 2.CePles2s€eTse3s | 
Dey 13.020l.0 =) | 
| Ghut*2) | 
J eReeS_ee2eye(3.0 = | 
1 # 2.0¢V*lYY - I 
| y.C))) I 
" I 
| RALOS = I 
§ 0. 75C47*FCRE I 


———_—n eee 


Fisure & 





1.0 


in computer time. 

Function RHLDS, in this case for elastic-perfectly 
Mersolec material, predicts strictly an elastic type failure. 
Thus, the full pressure reduction from REDPR (see chapter 


three) is employed directly to compensate for imoerfections 


Poo resiaval stresses. 
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ee 
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AXISYHMEDRIC (YIELD) FAILURE 


Wiewrca! Core of the optimization program is the 


ot 


exisymmetric failure mode, for the program designs its 
scantlings to fail by yield (see chavter one). Three 

meme rorrems in addition to FRANZ, PULOS, and REDPR are 
Mmacinuded in this BLOONS: Cr ol Ce and teas (eee 
figures 9, 10, and 11). 

The rather famous Vom sanaden and Gunther formulas G2 
and 924 (utilizing the maximum shell stress theory of 
fankine at the frame and midbay, resvectively) were used in 
design for many vears. Sedeouelee however, meny more solutions 
Poe appeared. With the advent of S-P stresses!8 , the 
stress an@lysis alone has improved in accuracy. ‘The manner 


in which the stresses are used to vredict collapse by 


axisymmetric yield varies greatly. Generally sveaking, the 
maximum strain enersy theory of lilses and Hencky vrovidcs 
the best manner of stress combination to vredict failures 
within tne various shell yield formulations. Tne point at 
which this is applied is also subject to discussion, 
Althoush it is generally agreed that the largest stresses 
Pevually occur at the frames, it is becoming evident The 


Pepa indicate the best oredictors use the midebay oOrea as 





Pie tHiviation point of yielding failures. The most 
irequently accurate analysis for axisymmetric yiélding 
moiyre 1s that due to ene. Although his formulation 
hes Nov been tested through SOMES. FOE of geometries 

and depths, tne tests that have becn made indicate his 
Solutions are at least as accurate as any others (to within 
1% of actual failure pressure in many cases), and mucn better 
than 92 or 924. In reference 14, Lunchik shows that very 
Pmecesciul correlations were obtained with tests of ring 
stiffened cylinders ranging from Xr (thinness Vat uoueo. le ve 
A =0.70. He recommended his formulation, however, only for 
"cylinders where scometries are in the range of axisymmetric 
feretoins", precisely the case in this progrem. Basically, 
Lunchik assumes a standard three-hinge felilure me nee 1ism, 
Doatuiatvine that the frame plastic hinvees fail first, ene 
predicting the pressure at which the mid-bayv hinge is 
complete. The basic difference between Luncnik's snalysis 
end others is his computation of tnlascic reserve strength 
His structure does not fail when some outer hull fiber at 
mid-hay has reached fielomscress. tt faits only afver this 


ity has vrogressed through the shell at tnat voint 


tte 


DIlASvic 
fear enoush to vrodiuce & hinse and precipitate failure. 


FUNCPION PiNCa The use of Function HEDPH (see chapter 


three) requires hoth an elestic end an inelestic collapse 
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duh, 


mressure for the barvicular failure mode under examination. 
mor tne Mier: icsperrectiy DilasciLe Macerials, there is no 
Winelastic failure pressurc" as such, since no strain- 
herdening is involved. Thus, the viold failure pressure 


Meperineo 2bove is substituted. For the solution of the 


elastic axisymmetric failure pressure, there remain two 


pO flies. One is the “exact" solution offered by the 
S-P stress analysis. When GAMA21.0 the Slaseic axisynmetric 


| 


his pressure i 


6) 


meae occurs. However, the solution for t 
Temmeup in & new "fF function" and requires a rather 
involved iteration. Because the REDPR method is approximate 


Mmemaay Case, &2n exact elastic axisymmetric value is not 


meeered, nus, the S-F solution was rejected in favor of 
Punchik'’s inelastic axisymmetric anelysisl, aUon ets analysis 
Moovecy Similar to Heynolds' analysis, for asymmetric . 
buckling (see chapter four), ioe beCan pele pelied i 
the strain hoerdenins (inelostic cease) using the secant and 
tangent moduli. By setting ETAN=0 and ESnC=l, tne solution 
breaks down to one for elastic axisymmetric buckling. This 
is the pressure computed by Function PLNCK, and the vrocess 
is taken directly from reference 13. 

Pleo bic. The “inelastic" pressure used in 


REDPR (and also the pressure to be reduced itself) is 


Computed in this subroutine. As expleined ebove, Lunchik's 
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yield analysis for elastic-perfectly plastic materials’ 

is used. S-P stresses are used in the formulation (see 
fieure 10). This, of course, meant an iterative process, 
and occasionally convergence problems were encountered. 

If indut geometries were satisfactory, convergence was 
Mecomplished in four or five cycles. Occasionally, the 
iterative program (TH“NS) calling ELNCk jumped outside the 
ranzse of convergent geometries in its search process. For 
such cases, an iteration limit of 10 was DUG lO mea. 
This meant that the oressure Goins packrto Vahila was not 
mierrely accurate, but good enough to continue with a search 
pattern to find @ convergent geonetry. An additional 
"safety valve" was built into Subroutine PULOS (see chavter 
t10) to prevent GAMA 21.0 and thus producing imaginary values. 
et GAITA 21.0, GAMA was fe equal to 0.95, and the shell 
thickness adjusted to achieve this. Thus, depending on 
entering geometries, occasionally shejl thickness itself is 
eajusted within ELNCKs. This is justifiable, in tnat the 
overall program is designing to a yield, not a buckling 
failure. Any buckling geometry, even if it could be 
converged uvon, would not be desired. Junchik's "yield 
pressure", PY (reference 16, equation 35) is the result of. 


the S-P stresses obtainec, and is the pressure at whicn 


ry? 


vVieldins bezins at the outer hull fiber, at mid-ney. His 





hh 


Plestic reserve strength ratio, PCLEI/PY (FCTR) is then 
PempuLred 2nd multinlied hy PY to obtain Lunchik's final 
Soliapse pressure, PCL!l. The eae certain of 
Lunchix's equations, mede somewhat confusing by SOF ea 
fn reference 14, is done in detail in Appendix EB. Before 
sending the collapse eee ce VOmnbiao. lt LS reduced { On 
residual stresses and manufacturing defects by EREDPH. 
mmenis polnvl (see figure 12), it is interestine to 
see, at least in one case, how the lunchik and Reynolds 
analyses compare with the log-log plot of hoop stress vs. 
the Windenburg-Trilling formula vresented in reference 8. 
Bomeenic particular hull diameter chosen, Reynolds’ pressures 
follow Windenburg's almost exactly. Lumnchik's vressures 
sriow hoop stress, at least in this case, to be rather 
eonservativ 
Bees ke Tahto Mihis subroutine uses BLNCK jin an 


iterative vrocess to converge on an exact shell thickness 


Which will fail by axisynmetric yield at the desired 
Pollacose pressure. if VHRANS cannot converge on a thickness, 


Mois Obvious that the input scantlings, despite chanses 


Mmimtiechness, are such chat failure by elastic axisymmetric - 


memewline occurs before axisymmetric vield failure. Since 


BOP Og Beate « og Chet 


Maem a Ssturensth of this metal is not then being used, 


PP. 


mee type Of Solution obviously is not desired. in such 
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Pecese, GunANS transfers control to tne next iteration loop 


of the main program (see chapter ten). 
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CHAPVaAR SIX 


HULL, PHICKNESS JUCREASE 


oa. ‘om er Ce Be Ble ated 


menerally socaking, the safety factor introduced by 
REDPR (see chapter three) to compensate for "imperfections 
and residual stresses" vill be adequate. However, it should 
be remembered that the collapse vressures develoved thus 
far all are "trigsered" by stressvalues at midbay. The 


Meee st stresses actue@lly encountered usually occur at the 


yer, 


memmeeteyine "flance", To account for this, stresses in that 
area are limited to 75 of yield stress at overating depth 
(assumed here to be 2/3 collapse devth) by subroutine 
OlRTHK, @s reconmended in reference Ewe iors Cone vO 
account for such things @s low-cycle fatigue, creep, stress 
Merrosion, and to insure & reasonable stress level in the 
Mmeeme Tlange prior to collavse for those frames with an 

; 10 
anivial ovteof-roundness  ,. 

Subroutine STRTHK (see fisure 13) uses the S-P stress 
analysis! ® EG COnpwle 421) Fe stresses of interest in the 
wiewiear tbhne fremes: inner and outer plate Beericee:. ont 
Honsitudinal and circumferential stress directions. If the 
largest of these is sreater than 75 vield stress at 2/3 


collaxse pressure, Shell tnickness is increased in 


increments of 5% until the criterion is met. 
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CHAPTER SEVEW 


GENBHAL INSTABILITY 


dne problem of findings accurete collavse pressures in 
the general instability mode has generélly siven analysts 
more Lrouble than the first tW0 failure modese The most 
accurate elastic peneral instability failure analysis was 
Mime by SoS. Kendrick in 1953 at Britain's taval Construction 


( 


KResearcn Establishment, and is eenerally known as the 
"Kendrick Part III" eoreone.. One year later, A. R. Bryant, 
working in the same ee iierment. developed a far simpler 
approach! (Hendrick's was exceedingly complex). Bryant's 
solution was a single two-term equation which could be (and 
hes been) used for design studies without extensive 
Pempuverization, although its accuracy left something to be 
desired (it was non-conservetive). wasically sveaking 
feyeane's equation incorporates the "split rigidities 
concept", where one term redresents the contribution of the 
Srell, and the other the contribution of the frames and 
frame-lenzth of shell plating. Although kendtrick's solution 
fs put in a simple gravhical form by Reynoldas?? (later 
extended by Beale jy, 1G wes rejected for use in this progren 
for two reasons 

1) since the general instability pressure is used 


merely to check the solution designed for yield 





failure, extreme accuracy was not reauired; and 
2) Kendrick's solution is good only for elastic 
fealure, And cannot be app*ied to strain- 
hardening materials. DI nCce@ol Sl DOr ran, 1s 
designed to be easily converted to one which 
can @lso handle strain-hardening materials, 
Kendrick Part III was rejected. 
eryent's solution was also mMejecvco TOLr the Gecond reason 
above. ilore recently, & very convenient and more accurate 
formula has been developed by Krenzke and i Aege- Geis 
very similar Mi Seructure ve Bryant's formula, but can be 
Mseqd for cither ideally viestic or strain-hardening materials. 
At about the same time, a very similar solution was worked 
out by unchik at the Model Basin. Both the Krenzke-Kiernan 
Pmaecne lunchik analyses give about the sane results when 
compared with actusl test aata®. Krenzke~hilernan WAS 
Pemecwnat erbitrarily sclected for this vrogram, solely 
because it appeared to be referenced more often in the 
imrver2ture. 
eee eCULIIN Baus Wis subroutine is called by tne 
meneral instability subprosren, and essentially computes 
failure pressure vsins the Krenzke-hiernan formula, It 
Blso returns the value of "li" (number of circumferential 


collapse lobes) which gives the lowest (nost Eruurcal) 
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mmure pressure. it might be noted that this program, like 
eeme OLhners in the optimizetion, is designed for dual usage: 
by ideally plastic or by strain~hardening materials (see 
fisure i/)). | 
SUBROUSINE GINST (see figure 15) GIUST tests the 
meeonins Scantling set from the main vorogram for failure by 
general instability. Depending on whether the collepse 
Peecsure is too shallow (test failed) or too deep (design 
too conservative), heavy frame (i.e., bulkhead) snacing is 
adjusted so as to give a collapse vressure that is either: 
1) ereater than 1605 times desired collapse depth 
(this requirement for a srall safety one is 
due to tne uncertainties of general instebility 
design, and comes from reference 10), or 
Poatess then 2 tines desired collapse depth. 
If, in this process, bulkhead spacing becomes longer than 
ie Chuire pressure hull itself, BS is set equal to the 
hull length (HULNTH), and the requirement for heavy frames 
is dropped (see chapter eight on Subroutine SVYFRH). 
Incoming bulkhead spacing may either be declered or 
undeclared in tne input data. If undeclered, starting 
bulkhead spacins is taken as approximately twice the null 
mieaneter, but always a@ multiple of tne small frame spacing. 


Tiewirede term In tne erenzve=-kiernan forrula contains, 
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@s one of its variables, tre “effective moment of inertia", 
Mme ints 1S cefined as the comhined moment of inertia of a 


Memeeeiceme 2nd its effective lenstn, EL, of snell plating. 


aiid Ot OE eee Sie Leen ied otal te ee 





Mais differs sligntly from Pecewemvene formula, wAicn mereiv 
uses a frame space (FS) of shell sil une hac ibs 
Getermined vie. the S-P stress analysist®, using the equation: 
EL=SLsPit3B 

jae use of eivher HL or FS in tne formula makes little 
feeererence in most cases; nonmenrere Sincessiae the tools wexe 
handy, the EL was used wnen possible. Problems of 
convergence, however, e@#ein develoovead in the PUILOS subvrorram. 
Wnen general instability voressure became very large (and it 
does for many geometries), the value of GAMA (see chavuter 
two) exceeded 1.0. Thus, a separate subprogram (PULCS1) 
was added, which set SL=FS when GANAS1.0. This is Oe lay 
acceptaole, since this case occurs only when the general 
Merete bprlity collapse pressure is far too deeo to worry abovt 


(its accuracy is only slightly degraded anyway). When 


td 


Bem@ecessive iterations of GINS? exvand BS and bring this 
collavse vressure-to shallower depths, 2L can again be 


accurately comouted by PULOS1. 





eee bes 8hGr | 


HEAVY FEAMES 


Subroutine HVYFSN (see figure 16) computes the minimun 
mzey Ot neavy frames needed to insure the general instability 
collapse pressure calculated by its preceeding subprogram 
GINST. Until. recentiy, heavy frames were designed using the 
ereandard Levy LOGIN ower es TOUNnG after considerable 
mesving, however, thet this formula often gave unsafe 
memetes |. A new formula was derived at the Model Basin 
by Blumenberg in 1965, which agrees much more closely with 
Pests. tesvlts from testing indicated that the effectiveness 
Memeo Licular size of heavy frames decreases as the . 


cylinder is lengthened and @lso that the minimum size of 


heavy frames necded to locelize failure between heavy frames 


° 
e 


is possibly not dependent upon their spacing”. Although 
adequate testing has not yet been publisned te positively 
Pemrirm Slumenberg's formula, initial results show it is 
better than wnat was formerly used, and thus it was put into 


ey Yi: 


7 
ne et -, where: 
CChve2-1.0)eETAN) * 
Wmeevey be subeeituted for ETA with idesliy plastic 


Peer tolews fee Sctualily is not in Blumenberg's 


formula, put was vlaced there by the autnor to 


Peso iy. a eer in an inelastic cptamization. 
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= - 2) M is the critical buckling mode for elastic general 
instability failure of the cylinder with the neavy 
frames replaced by typical franes (i06., & cylinder 
equal in length, but with no heavy frames}. 
oe cowec! soyim5ols” for other variables) 
imimecompitings PC, the radius to the combined centroid of the 
heavy frame and its effective length of shell Pees an 
ELH had to be determined. This wes also vrovided in 
Peumenbers's report: 


a ee SLs FIs (AF+S -SL* ) HS 18 
~ APR+SLET 


1G can be seen to be a form of simple Poconos O01 one 
original EL formula developed in the SP ete. 

Another feature of HVYFEM is to send a value of zero back to 
the main program if the heavy frame snacing computed by 
GINST is equal to (or ee Eien ethe enn) lencth.”. This, 


of course, would mean that heavy frames are not required, and 


Theat the hull end closures a&re used in their place, 





ChAPTER NINE 


oom =D SrA Ce nbn) Hine 1 


A 


The simple weight/displacement ratio was selected for 
Paeroptimization oriterion. Not only wes it the ceimplest 
moeuse, but it also seemed to be the most general, all- 
encompassing determinant of On optimum design. Due to the 
subprogram system utilized in the optimization, another 
form Of Criterion could easily be substituted if the need 
BOSC. 

Essentially (see figure 17), WEDSP computes the weight 
of a hull section, the lensth of which is ecual to one 
heavy frame spacing, splitting the heavy frame on each end 
in half. This weight is divided by the same ‘pon S 
displacement in sea water, taking into account internal or 


extern2l frames. 
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CHAPTER VEN 


HAT PROGRAN 
Esse ney the main ovrogrem verforms very few 


. 


@eculecions. Its primary task is manipulation of tne 


a ad 


Verious subprograms and iteration control of the entire 
Meocess SO 42S to @rrive at an optimum design. See figure 23 
for & basic, easier~-to-follow main flow diesranm. Figure 24 
is the conpviete main prosiam flow diagran. ech en OM, » 0G 
computes "reasonable" entering scantlings for the iteration 
process, based on input data and proven desisn practice. 
dne objective of this tyne of main program design wes to 
oeain 2. system ijn which the various modes of failure and 
safety factors could be chansed or interchanged easily 4s 
desired. As all enue modes are separated into 
Subprograms, it is possible even to substitute an entirely 
Garierent enalysis for any individuel failure mode. safety 


mMecvors, which so often are subject to modification due to 


> 


mew test date, are also easily revlaced or cnanmed,. 

[he input data required (see Aprendix D for formar ) 
are: overall hull length, bulkhead spacings (ray or WAY mom 
Mema@iven. if not, 1+ computes on optimum BS), internal or 
external framing lesired, founs's modulus, mecel Hail Secoss, 


a 


Meeeeedicisitv, 2nd Poisson's ratio. In addition, as presently 


7 


Berouo, tne program mey be iterated for verious innut values 





of depth and/or hull diameter. 
The depth input (in feet) is converted to a collavse 
pressure PC (in psi) requirement by using theequation of the 


OO. 
? 


mean line drawn in fisvre 1 Chis represents On averare 


@epour vs. psi curve for thes verious oceans included. ‘he 
figure (minus the mean line) was taken from the handbook of 
Seeen Hneineering ables, published by the U. S. Naval 
Ceeanogsraphie Office, and compiled by £. Tee eke Tere 

Pie ixrame constant CC is computed from the equation of 
the curve shown in figure 19. CC is used to determine frame 


proportions to be used for various depths. The curve is an 


Sverare Of proportions of many ring stiffeners used in 


present generation subnersibles. ‘The method of ohtaining 
frame proportions is printed in the program outout (see | 
Awpendix D). <A more detailed explanation may be found in 
Bppendix A. 

mie prosram nexvy computes its "detum voint", or 
Meemenic Set of hull scantiliness. Shell thickness, '', is 
computed from the simple hoop ye formule: 


Oto eC. <DEL 


@7~ 202m & 


an O#SiGy 
PreamenSoaAclie ys Or, More accurately in this case, 


unsuvosoorted lens ta of shell between frames (SI.) is computed 


X 


from the varemeter §, where: i evanersy 
: oe AC] =) 


oi 
me 
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Reference 6 Lists the usval range of 6 bem present doy 
Beupmersibles as 1.0-2.5. ‘The main program (see below) is 
set up so that it starts at a maxmium freme spacing and 
then decreases it in further eee Jee wan ce) Boieal 
Seerving value of Q would be desirable. After several 
calculations at various devths and geometries, a value of 8 


? Ad 


moo was vsed for the datum point. This gives 2 wide range 
of SL valves. The iterative process of Ais pea: produces 
values of 8 that eventually go low enovgn to bracket the 

eoove renze, Other methods of obtaining a starving SL were 


) 


invesvisated, such as the combined solution of hoop stress 

end the Windenberg-Trilling equation, but all gave too wide a 

france (usvelly too large a frame spacing) as depth increased. 

Thus § ; which does not vary greatly, was selected, 
festarvinge frane size is selected based on the ratio of- 

frame cross-sectional area to the cross-sectional area of 

one unsuyported length of shell pleating. Reference 6 lists 

the normal range of this ratio to be 0.,2+-0.8. During 

rather extensive investigation, however, it was found that a 

meacio of greater than 0.5 gave frames that were grossly 


Overcdesigned. ‘hus, the starting vatio vas set at 0.5: 


~ 


Ar. teal s| “efit 
Sud gue A 


= 
ps 
i 
© 
WW 


ine frame moment of inertia is then computed in accordance 


merci relavions develoved in Appendix A. Criginelly, it 





vi 


wes intended to select a starting frame size by @ more 
accurate approacn. However, the only formula which could 
Bewsolved in anything approacning a closed form was far tco 
Bonservative (i.e., Tokugaua's Cone, reference 25). The 
eadventage of & computerized approacn, i.e., investigation 
was the deciding factor in using the more random iterative 
method described above. 

mt bonis point the program goes into a double loon, 
Mmeeracins on SL (outer loop) and AF (inner loop). ‘The 
Seeiulings are tested by HivLDS,. If they do not faii in the 
lobar buckling mode at design collapse depth, the vrogrram 
Pees On to THanko. If the scantlings fail BeLDS, SL is 
decreased by nultiplying it by 0.9, @nd the RiuLDS peat 1S 
[erui. the decision to run RBiLDS at thls point was made in 
order to start as near the "shoulder" (or at least, not to 
fie right of it) of figure 12 as possible, since this is 
generally acknowledged to be an ares. Of cCpUinum design. 


ry 


mimeo, the decision to vary SL, rather than 7’, to achieve & 
Meio, Scantlingss that would not fail the RNLDS test, was 
Beage for two reasons: 
1) From figure 20, it can be seen that either 
Pomme mecOulLGmbo SUCCeSSIUuLtvensed cto change 
nn 


Pies atlure Dreceure, vivh 1 being 


slightly more efrective. 
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2) The next subprogram in line is THKNS. it is 
sonnel Be ANLDS vie enother iteration loop. 
Yo have both programs converging on different 
valves of T would cause endless loops in almost 
every case. Thus, RULDS was iterated on SL. 
The decision nov to vary AF to change RXLDS feilvre vressure 
Hs justified by figure 21. The Nott buckling equation is 
almost identical to ENLDS, and thus can be considexed 
the same for these tye investisations. 10 may be 
Mmeved that for some PeOmeunies: me has meowefiect at 4ll on 
PCR; thus, AF was not used in the RiLDS loop. This sane 
reasoning (see figure 22) vas used in deciding against using 
Pome Ciat ion Of any Kind for convergencds within vhe 
TIKNS SUDOrOcream. 
As stated above, once the scantling set (revised or 
mareyvisea) gets successfully past RNLDS, it goes on to 
PHENS. There, the hull thickness, T, is adjusted so as to 


heave the shell feil by axisymmetric yield exactly tenia 12) 


a 


at design collapse vressure. At that point, the scantli 
Min revisea tT are again looped back througn HNLDS to insure 
against lobar buckling, and SL adjusted again os necessary. 
ieeol has to be re-edjusted, the scantlings egain go through 
Paniico. The vrocess continues until either: 

imouto UNOl Myre leme Sceyncgilanets vass without 


cnanse throuszn both RALDS and TERNS or 





COLLAPSE PRESSURE , PSI 


GOO 


600 


£00 


EFFECT OF FRAME SIZE ON ASYMIMETE 


BUCKLING PiREOSU 


2 


i om 


Sy = 88 000 pst 


BYOtie,. 


Ge 
21C 





0.2 O*+ 


& —_ 





an tae Eo 


0.6 0.2 


FRAME AREA/SHELL AREA 


(NOTE: TAKEN. FROM REFERENCE 15) 
une gan aly iat & 21. 


10 


ET EP EE AY 


MGINLEY 
4-770 





FAILURE PRESSURE (PS)) 


2000 


1000 





EFFECT OF FRAME SIZE ON AXISYMMETRIC © 
(YIELD) FAILURE PRESSUPE BY LUNCHIK. . 


MAW Lolly SO 
HULL DIAMETER : 16 FT 





DES OS ea Ee. Se ee PS Te. 


5.0 10.0 15.0. 
FRAME AREA CIN7) 


Pecure 22 


MCGINLEY 
4/70 





2) SL becomes so snell thet there is less then 


Pour inches clearance between adjacent freme 
flanges. 
On tne iatter case, the loov is skipped wi shout printout, 


— 


ana the next iteration is Feared, Much the Same as what 


ad 


@eemrs tf LHKNS cannot internally converge on a shell 


thickness (see chapter six). 


t 
Je es is Loe | 


Once the scantling set gets vast THHNS, VY is again 
Pojtoaved Wovard, if necessary, by STRTHK, 4. loop to re-test 
back through RNLDS and THKENS is not used here, and would 
Semverno useful purpose, Since V's from THKNS and STRTHS 


could rarely ever be made to converge (l.e., the lergest T 


meee 
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from DOun oOrogranmns 1 

Prior to entering GINST, the integer L isset equel to 
O. This, used @s one of the entering arguments for GINST, 
ms 2 Control variable. =O insures normal (iterative Ps 
fesie@n) operation of GINST. T=1 is usely solely to obtain 
@ single-pass value of i] for HVYFRM (see chavter eight). 

The scantiing set then enters and is tested by GINS, 
io 1s Aaadjusted @s necessary to insure £©CG falls between 


= 


Meo fro end 2PC; the only exception being if DS becomes 
sreater than total hull lenzth (see chapter seven). ‘The 
meeeow tna So, racvner then 1, SL, or AF, is used here to 


Ponverze on & Satisfactory PCG, is because in virtually any 





ai 


‘ease, regardless of depth or geometrical vroportions, BS 

has a much smaller effect on hull weight for the same PCG 
reduction or increé@se. Since hull weight already had been 
optimized for shell yicld failure (the desired situation), 
any further wv, SL,or AF changes would most probably take the 
design seriously off the optimum. 

Once past GINST, the secantlings pass through HVYFRI: 
to obtein a suitable heavy frame desisn end WYDSP to 
emovlate that varticular scantling set's relative efficiency 
ee Oz. 

AG this point, values of the scantiings are vrinted out. 
Generally speaking, the print-out of each Line of scantilings 
can be said to be optimized on thickness of the shell, 
elthough certainly other factors (SL, BS, etc.) change as 
necessary to keep the design on & yield failure besis 

Cnce out of the first loop cycle, AF is decreased by 
maicipiying by 0.8. This is done for ten cycles, so that 
the original AF is reduced to 0,108"4AF in the lest cycle. 

SL is reduced in the same menner in tre outer oe and @s 
described above, also may be reduced within eacn iteration 
as necessary to vass Ri LDS. imethneeecnd os boun oops (100 
major iterations) the program terminates for that set of inout 


Parameters, &nd prints out the optimum (i.e., smallest 


Weight/disvlacenent ratio. 
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CHAPTER HLEVIsn 


Bis HOU es) Als ‘Dp _CONCLUSI Oho 


Various parametric studies were performed with the 
orogram, once jt was checked against some contemporary | 


Piomersi bie hulls to see if it indeed was "in the ball 


Sy) 


park", The most obvious study was to see how the W/D 

ratio varied with depth (figure 25). This was done for 
three common steels in use today, &@l1l with PN eEaers Gey Sala 
Mideally plastic" stress-strain cheracteri er The points 
obtained plotted into smooth curves on the seni-log: plot 
meeo., tit is fairly el that if a W/D greater than 0.5 
is considered unsatisfactory, the following limits would 


BDIOLY : 


Suet pune OR N/a 5.0 
HTS 3500 FD. 
HY-80 6200 FT, 
HY~100 7400 IP. 


A second study was conducted to observe how the WAD 
Peervo variec with hull diameter, @ll other factors 


Miscludine depth) being constant. it was found that it 


— 


Peerca very Jittle if eb @l1, as can be seen from the 


POlLlLowine set of dates. 





WEIGHT/ DISPLACEMENT RATIO 


OPTIMUM WTI/A RATIOS FOR 
VARYING STEELS AND DEPTHS 






CONSTANTS: 






DIAMEstein 05 SOF I 


HULL LENGTH: 200 FT 





FRAMES * INTERNAL 





\ (2 5 4- > 10 


DEP UrCUSANDS OF FEET 
Lees = 
MoOGINILEY 
é-/10 
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Given: HY 80 steel, collavse denth=5000 ft.: 


HULL DIAMEPER (FT. ) W/D PATIO 
AAO, O.417 
30 0.418 
50 Oren 7 
80 0.409 


It could tentatively be said then, that for a given depth, 
hull efficiency (W/D) will be approximately the same 
rezardless of how large the hull diemeter may be, 

Figure 26 shows the result of plotting the variance of 
feewerien metal yield stress. This also plotted into a 
meevr Curve in this case, i1t is probabie, however, that 
for higher yield strength (and therefore, strein-hardening 
foes eCUuLTing &2 slightly different analysis, the curve . 
would have a@ sharp break. 

iieecnird plot @ttempted (figure 27) at first appeared 
to be & hopeless scattering of data, bul atvcr Sone anatysis 
revealed rather interesting results. The points plotted 
were taken fron & single optimization (one diameter, one 
depth) such as in the example vrint~-out in Ax ppendix D, A 
CBliige computing Uhiinesomfacvor Was inserved in the progran 
and vrinted out with the regular data, 

The "SL loop" is actually a series of W/D points 


Bemonted Tor scenvline Sets with the same Frene sracine, 





WEIGHT / DISPLACEMENT RATIO 





Be 
OPTIMUM WT/A RATIO VS. STEEL STRENGTH ~ 


CONSTANTS : 
COLLAPSE DEPTH: IGOO FT 
Fier eee PAE iE Re Sie) Sih 


HULL LENGTH: 200FT 
FRAMES: INTERNAL 





40 60 80 JOO 


YIELD STRESS, THOUSANDS OF PSI 


Fissure 26 
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Puy verying frame size as an input. In each case, it is 
strongly evident that a saddle point, or optinun W/D for 
that frame spacing, was reached. ‘This means That for at 
meast this particuler set of inputs, the program's frame 
size (AF) iteration ranze was large enough to bracket 
optimum values. One of these frame Spo UME ninimums, then, 
was the optimum W/D ratio. 

The "W/D envelope" encloses all W/D values computed in 
the program. There seems to be good indication that the 
Nootimum W/D" indicated is a true ovtimum, since the lower 
DOrvion of the envelove rises on either side of it. 

Mother indication siven by the plot is that the 
program gives a much wider variation in W/D ratios with 
JjJargzer frame spacings, for varying frame sizes. As frame 
spacing becomes smaller, the W/D ratios produced become more 
lieonversent". 

Perhaps the most obvious conclusion from figure 26 is 
that A eWionecuwmemCerusiniy NOL an accurate indicator of 
optimun W/D ratio, although it could be utilized (after 
Petenorve data gathering) to indicate the general erea in 
miowenm CO aesign. 

it ts interesting to note that average computation - 
time for e@ch optimum computed was less than one-half 


Paves GSO, 20 Was Anverestins to note trat tne 


— 


Tl 


MS citvin soluetons contained @s scantlingzs, generally 


88 


speaking, smaller frames, larser heavy frame svacing, end 
Slightiy thicker shells than submersibles of present design. 
Evidently, this results in better W/D ratios. 

One of the major advantages of this program and its 
Benerel method of computation end print-out (see Anpendix C) 
Mmemoiet it Sives & great variety of alternatives from which 
to choose. .For instance, perhaps the optimun W/D ratio 
Remearcervain hull configuration and dceoth contained 


& 


scantlings which gave a very small frane spacing. 

may oicolly, this may give a superior W/D retio. But 
Peomerca) ly speaking, the cost of fabricating and installing 
& great number of frames may be out of the question. Thus 
(pearcicularly if the submersible is not critically weight- 
twinced), the table of printouts may be “browsed" for 

more attractive scantling sets: ones with acpawilinle W/D 
mepnosg, but with inherently lower construction costs. lie 
Should be repeated here that each of these printed lines 


"| 


ere not mere rendom cnoices. Each line, vrior to pvrint-out, 


nas already been through one of the main hurdles of the 


= : 
ia 


Bobwization proferem. The design of shell thickness and 


frame spacing combination to give failure by axisymmevric 
yield, and thus most efficient use of the material's 
Mine reis SUrenstni, LS completcd prior to printing each line 


Sette ensiwer Tanlc. Thus; the voroszranm ray be used not only 


Peeeorrucvural optinization, but elso is @n indispensible 


reference for @ny economic analysis of a particular design. 
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CHAPTER TWELVE 


RECOM ARDATTONS 
One wrest iw A ee Oar te ONG Te eee =e 


a: 


fhis Program, 
1, Ranses of Variation. One of the most common 


GU sSulOms set LOr 2) Dorvicular run completion was 

"i wonder if this is meow bye Ulic Opulmumn, OF 

didn't I go far enough in frame (or frame spacins) 
variation?" lt is believed that plots such as 
feeure 2/7 tor each run could probably give a 
aGefinitive answer in most cases, Not only can 

Cue yvcll Mf cech See of frames ran throuch an 
optimum, but if should also usually be vossible 

to tell if the 'W/D envelove" passed throurh its 
optimum. This may be @ rather peinful way to 
assure one's self that his run covered all the 
Betti tory ia Geres necessary, but 2G this writing, 
imum PCAs evembece cic SUTESUL Way. As an alternative, 
some rather extensive studies could be conducted 
vsing various depth and otvner inouvts to test range 
validities. Once determined, the program's range 
of AF and SL variation, governed by the indices kK 
and J, respectively, could be controlled 
appropriately. From a rather cursory inspection by 


4.) 999 4-e. 5k. es psrieneor) lyn 4-4 ae So eee 
Cie w eer yp tu popearca bhavy 19 BUNS conductca 





Nw 
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in data-gatherine for this paper, the ranges of 
=K=10 were adequate in all but a few questionable 
igrecus Olugraurication Procedures. It is rather 
Sovvouowonal Wiviour any GCouby, the greatest 
Siler ecomude mig ENC ODGimMI1ZaUC10n 1S in its 
allowances for such relative "unknowns" as residual 
welding stresses, low-cvcle fatizue, shell or 
frame out-of-roundness, etc. ‘The gross 
approximations made by REDPR end STRTHK certainly 
fall far short of the accuracy of the various 
failure mode analyses. There is no doubt thet more 
basic research must be accomplisned in this region 
before & completely dependable optimization 
program could be achieved. As it stands now, the 
Sapevyeuercvors PULL anto. tne proerain coulda fall 
H1vO vO CalLecor les: 

a) Completely safe design, in which 611 
scantlings are overdesigned to the extent 
that desisn "optimization" is almost useless. 

b) more re@listic (i+e., lower) safety factors 
based on scanty exverimental date, which 
toe notemniverss ily eoplicable, and thus 


a ~% PAs wenn 
mieht be considered unsefs. 
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Jt is hoped that the safety devices employed in 
this vrogram adhere to a "middle of the road" 
Dee ee bes believed, results in 
somewhat of an overdecsign for the two between- 
frome Peviure modes, but in underdesizgn for the 
general instability mode. This underdesign is 
hopefully picked up in STHTHK, which, as mentioned 
in Chapter 6, insures "a reasonable stress level 
in the frame flange prior to collapse for those 
Reames mith sedan tia | Biemeoindnesst! 0, 

‘ne ONLY Way, it appears, to resolve these 
questions, eee extensive testing of models with 
deliberate, measurable defects of all types. The 
most needed data of this type is in the failure 
mode most effected vy defects: general instebility. 
It is believed by some euthors~ tiavecne Cul of 
roundness analyses developed thus far are overly 
vessimisitic when applied to full sized submersibles 

Future Programs. 
The mosc Obvious extension of this program is into 
strain~hardening materials. Because of the 
program's meneral: characteristics (i.e., @ main 
flow control vrogram manipulating subprograns 


which ectuelly do pressure celeulations), it can 





be easily transformed with.very few changes other 
Gislh Sibprocram additions. 

The need for such a transformation is obvious. 
when one observes hull steels projected for future 
(and in many cases, voresent) usage. Any steel 
with @ yield strength greater than 100,000 psi 
can be considered to be @ strain-hardening material. 
Detaiis of the method of transforming the present 
program into one suitable for use with high-strength 


steels may be found in Apovendix E,. 
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APPENDIX A 


FRAME DINZHSIONS 








De, 


APPENDIX A 
Pee Dario t ONS 
Perea, FRAMES: 

iimordecr vO arrive av reasonable frame cross-sections, 
averages of a great variety of frame scantlings used in a 
number of submersibles at various depths were obtained (see 
figure 19, chapter 10). The common denominator for all of 
these freames is the web thickness, B. “his method was taken 
from Adamchak in reference 1, who did the same thing for 
surface ship frames. 

The ratios used in figure 19 were taken from the 
following, derived from the averages computed (see figure 
28): 

Flange width=Fw=CC+3}, where CC=y of fisure 19 
Flange depth=FD=PW=CC#B 
Flange thickness=1, 7B 

In eddition to the above, the remaining terms in 


foeure 26 ere defined as follows: 


ine Diciomeecmonehreine Cencroln trom plate neutral 
axis 
DN : Distance o7 combined centroid from plate 


neutral axis 
The various formulas for frame and combined frame-plate 


7 i 


tela thas c ua ef - on tL 5 we ie 2 fe aC ye ._2: = oe , ais 
Mertens Of inerti> ane for cenver of &Yyration were v4Nen 
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T- SECTION 
DIMENSION IDENTIFICATION 


DN 


| X— 
TS 
ails a 


Figure 28 


NCGINLEY 
4£(10 
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from 2 general development of these in reference 19. They 
are listed in easily prosrammed form, and according to 
mererence 19, differences from values of actual standard 
Broavetion frames due to fillets, etc., are always less than 
wWO per cent. 

The combination of these two methods cnables the 
computer to easily arrive at any variety of T-stiffenex 
meee veriStics, each of which is proportioned according to 
Piesenme submersible design practice. 

Peter SuvGtitucing vLhe expressions for FD, FW, TF and p 
into the equations for T-frames in reference 19, the 
following expressions result for typical (i.e., not heavy) 
frames, and are used in the program: 

AP=FUSTPTPDs5B = 2, 72 CCe Bue? 
Frome monent of inertia, Bis 
FJ = CCOPsBect, where: 
CCOP= (2,892 CCxxtt+11 , 34 CCue 344 3.47 CCe a2 
eee OC) (1 25.0-CC+20, 4) 


Effective frame~plate moment of inertia (using 


ah 
wenesee & SEO EPONA LP OS? FE LOT PGE 


errective Jenathn of shell plating, EL, develoved 


eee 6 EY Cet OR See ot Ey ye F—28_ 


inmemerence 18) ..HT: 
EL= 0,225%CCxBaaclbse (CCOHL. 7 )2242%((1.07%CC4+0. 562) / 
(OCH. 7)4+3.0%(1.634(2-1.7% Wey eo le 7) 40) eee / 


Ge ee rc se) /( that) ) 





HEAVY FRAHES: 
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Heavy fremes were averesed in the s&me manner as 


above, althoush less data exists. The relationshinvs 


obtained did not depend upon depth, and were determined to 


be as follows: 
PDH=17. On-BE 
eal Gis Orrell 
CRH=A, Ox Bi 


CCOP=1 330.0 


ty 


leavy frame area, AFH: 


POPPE, dee FOP en Oe me Oh gg Pe OE RR & rae EE Ot 





APH=443, Ox BH a2 

Heavy frame moment of inertia, PIMs 
FIH=1 330. OxBerclt 
PCH=144525.3H+1/2 0 


DNH=FCHsAPH/ (APHHELHeT), where 


ELH=effective length of shell plate next to 


heavy frame (see chanter 8) 


Bli was not needed in the pyrortrain, 


che 








APPENDIX B 


LUNCHIK'S PLASTIC HINGE ANALYSIS 


a7 





IOC 
APPENDIX B 
DUNC pio sO eC HIG AMALY SIS 
(Refer to reference 14 for this discussion) 

Lunchik's final Pc/Py equation is developed through the 
use of his parameters Band Ek. There are not solvable except 
by assumptions which Lunchik makes. Two of his assumptions 
ee ee; konuk, end Oma a K./Kg) are easily worked ovt from 
substitutions in identities from his paper. However, one is 
tess clear (i.e., Ovghy oF SOBs), particularly because Bg 
is misprinted as 6g in reference 1. 

Prom Lunchik's pee ons. DiiesmeneivelS Cf 2 One 
vnit square element: 

(1)Mekgp = circunferential edse moment 
(2 )Ng= Kyph = = circumferential compressive membrane 
1orce 
p Ke — 
(3) Bg= h Ke 
(4) Oz 


pe 


Kp = circumferential membrane stress 
Assuming (1i.e., approximating) the circunferential 


bendinz stress to be elastic: 


OY 
= 


_ 2e ia _ 6Mg 
Ouae “T= Ep Ms 2 LE a L2 (5) 








Fron (2) and (4+): 
Ne 
Nz Sail =) Cine si fe (6) 
Using (5) and (6): 


Sig SMe sh. Slop. Ske. 
Ovag ly No h: Koh hi Ky 





>) 
S 
TO 
! 
0) 
“ 
0) 
WwW 
a 


+.) ra 
“Eo 44 6 
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Hemiurvuier re-agtanwe these @ssumptions to obtain 
expressions for B, and By to use in vrogram ELNC 


O 
(7) Bg = s YS... . (8) K./ Ky a bun fe, 


Sie hK,B RK, (6 °%,-) 
< 


| Orgs “a 
_Kx ko = lr K (S Koy) ne bd 
BL ik, RK aK ii Fae) (9) 


(7), (8), and (9) are then put into Lunchik's equations 
12) ; (13), ree (hj) ous 0 O, and 0, , and these in turn 
are substituted in his (25) and (26), which are used to 
solve for the plastic reserve stength ratio Pe/Py 


(FCTR in progran ELNCK), 
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APPENDIA E 


METHOD OF PROGRAM CONVERSION 


TO ACCOMMODATE SOTATNL-HARDENING HATERITALS 








1.30 


Ae lex 
CULTDE eOr =e wnOD PCH COMVE ee 1 O aryin 


OPTIHIZATION PROGRAM FOR STRAIN-HARDENING HAVERIALS 


i 


Strain-hardening metals differ from clastic perfectly 
Meee Mevals, in that above the yiclc point, the stress is 
not a constant value as strain increfAses, but increases 
(usually at progressively slower rates) as stra incre®ses. 
| This means that true plastic flow is never achieved; the 
metal continues to retain @ modulus of some value, albeit 
sm@2jiler than the original constant value. This results ina 
type of combination elastic-plastic buckling feilure in 
peeeeeoned Cylindrical shetis, termed inclastic failure. 
fee particularly important to analyze subersibles 


A. 


constructed of strain-hardening naterials Migs Wile Sore 


BE 91 20. liormally, a submarine hull is 


hardening analyses 
designed to have some plastic yielding (1.e., initially 
beginning in the hull adjacent to the frame flanges) 
soimewhere between opereting and collapse depth; indeed, there 
must be yielding prior to collapse depth in order for the 
mimi to crush theréewe With ideally plastic materials, once 
mMiesvyiclains polnv 3S reached, generally sspeaking, buckling 
foeore 1s ruled out Cor 1t would have occurred earlier, due 
to hull geometry). With strain~hardening materia 


~ 


moverm@:, Ccnice ti® luli bosins to plastically defors 
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Meewtins tallure is not ruled out. In fact, there is a 
strong possibility that the buckling feilure (which will 
occur) will havpen at a lower pressure than the yield failure 
caleulated by an ideally plastic analysis, due to the 
Meemeei neval modulus. Thus, to use an ideally plastic 
collapse pressure analysis on & strain herdening material 
fede SLive Cangserously overoptimistic failure predictions, 
Peeerculariy if the metal has a very high yield point. 

For convenience in analyses, strain hardening materials! 
stress-strain curves are characterized not only by EF, but 
@lso by Ey (WTAN, tangent modulus) and Eg, (ESEC, secant 
modulus). See figure 32. 

Because with tnis type of stress-strain curve, the 
moduli are always dependant upon the stress state when 
meeve Vield stress, it is necessary for all strain-~-hardening 
collapse pressures to be computed using an iterative process. 
it is generally the approach to the solution of this process, 
and an example pressure analysis by Reynolds~° whew ows id 
comprise the rest of Appendix KE. 

Pesenvcie liv, che crivical collapse pressure IS (Ol onysliLele: 
when the buckling equation (depending upon buckling mode, 
references 31, 13, or 20) is solved simultaneously with the 
pre-buckling equation (stress intensity as a function of 


3 Peahere eA A iii 
oressure), where: 


Figure 32 


TYPICAL INELASTIC 






STRESS-STRAIN CURVE SHOWING 


EODULI VOR SINTER OT 


rreure 33 


INELASTIC BUCKLING 
EQUATION | 


PRE-BUCK LING 
<-— EQUILIBRIUM EQUATION 


eee COMLAPSH PRESSURE 


GRAPHICAL DETERMINATION OF INELASTIC 


BUGCERILIWG PRESSURE 
KSGINLEY 
Alo 
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Stress intensity= oF aCe =AOUC: See figure 33) 


Bither V-G or S=P stress theories could be used to calculate 
Mm@ewctrcss intensity, S. V-G is less accurate, but 5-P 
might give convergence problems, depending uvon hu2 
gecometry and depth at which stress is calculate 

One difficulty, that of finding 4 way of describing a 
strain-~hardening stress-strain curve with a minimum of input 
Gece, iS solved in reference 23. in this method, the entire 
curve may be approximated with extreme accurecy by using only 
four inputs (see figure 34): E, &. Oe and i, Shy 
manipulatins some of the Romborgz-OQsgzood equations, it is 
possible to obtain Ey and Hyg, given any value of stress, via 


the following expressions: 


n= jo OR NOS) 
1084. 0(O. /Os) 


B/E,= 1+0.12857nG/o,)”" 
B/E= 1+0.42857(6/g,)™ 

This will be @ssumed to be the content of a subprogram 
called Subroutine ZOHOS (see figure 35). 

Both the pre-buckling equilibrium equation and the 
buckling equation willl be approximated ite sen recon CF 
Mperesu With gurssecme dines. nis, and some of the followir 
methods of determining the intersection of the two equetions, 


a“ 


were patterned generally after similar metnods used in 


0 


rmererence 2 


1 34 


METHOD OF OBTAINING ROHBERG-~OSGOOD 


INPUT PARAMETERS FOR COMPUTING Eg AND Et 





Figure 34 
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ROHOS 
CoG cAN . 
ESEC) 


COMMON /G/ SIGA,SIGB 


EN=1.0+0. 3853/AL0G10 (SIGA/SIG3) 





BNI=EN~1.0 














BTANSE/(1.0+0.42857uEN« (SIG/SIGA)##EN1) 





























ESEC=E/(1.0+0.42857+(SIG/SIGA )~sENt1 ) 





PLOW DIAGRAN: 


SpUbEeU MWe ROnOs 


Figure 35 


MGitley 
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[ade al 


subroutine LINC. This subroutine computes EP 


PORES Bort Cae © iA Gt IP <2 gate Oe OPPO EaD 


line slope, and BEZ, the line stress intercept of the pre-- 
buckling equation, using the S~P analysis (see figure 36). 
The values given by LIVE will change only with scantling 
‘chanses. This subprogram would censtitute the only deviation 
of the main flow program from the perfectly plastic case. 
Jt would be placed in the main vrogram directly after point 3 
Seeeey brior to the HHLDS calljs It would also be used 
within VHRLINS whenever T changes. 
. MiecwrOllowing discevesron will anvoive the solution of 
the asymmetric inelastic buckling mode @s developed by 
Reynolds in reference 2Q. The same generel iterative process 
would be utilized in the solution of the other two modes of 
Moe fariloure. 

fipee vi Om PCH Secet ieure.3/, This sme 11 subprogran 
merely computes asymmetric inelastic buckling pressure, 
meine as inouts HTAN and ESEC (computed by ROMOS), T, SL, 
and PCRE1L (computed as PCRE by the elastic portion of FPCRP's 
calling program, PNLDS). ‘The pressure 1s computed using 
Reynolds! equations outlined in reference 20 (for 2 
simplifiec presentation, however , see reference 12, which 
gives the axisymmetric mode also). 

PUDrole een Gece Ticgure 93, This subprogram obtains 


Gnhe intersection of the pre+-puckling equilibrium equation 


> 1 ’ - » <a AGP be aa acy se : eyes pea : { ae ie 
Per eoe iene erica OF tne two input Dressures \ruisk 


(SL,T,FI', 137 


COMMON/D/E ,GNU,DM,4,CC,SIGY ,CCOP,RHO/H/BEE, EMM 





_/U/THETA ,ALPA , BETA 
P=0 e 9xPC 
CALL 


FRAME 
CFE a vol) 


SIGU=-PsRM/T 











Chic 

PULOS 
CDE SL Sas ag See Sa 
Pe Oe eae 


SIGNP=SIGU+s 





5xS1GU 


S(I)=SQRT (SIGNPs++2+SIGNHXs22-SIGMPsSIGHX ) 






ENM=(S(2)-S(1))/(P(2)-P(1)) 





Pero 4) laa ( 1) 


WiLOW seas SUBsOUTING LINE 


Figure 36 
MACCINL EY 
4/To 










POR 
(PCR Eee a 
ESEC ,ETAN, 
SL) 






COMMON /D/E,GNU,DM,Z,CC,SIGY, 
COOP uo 





NUP=0, 5~ (0. 5-GNU)#ESEC/E 


PCRP=(PCREL#(1.0-GNU##2 ) “ETAN/ (1. 0-GNUP«#2 ) ) 
(1.0+(0.922+SQRT(RN«T) /SL)«(BSEC/ETAN-1.0)) 


FLOW DIAGRAK: 


RC LON PCP 
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. 
RNPT noe 
(PCREL,T, 
Siu, Phsi.SiG1. 
Slee tee See 
P,S) 


COMHON/D/E,GNU,DM,Z,CC,SIGY ,CCOP,RBHO 
G/SIGA ,SIGB/H/BHE , EFM 


CALL 
ROHOS 
(SIG2,ETAN, 





ESEC) 


RioZ=eOhigGhe denon C man, ol) 





ENN=(SIG2-SIGi ) /(PRS2-PRS1i 









ERE=SIGL-ENN*PRS1 


P= (CEE~BEE) / (SMH~-ENN ) 





S2:HNN+P+ChE 


FLOW DIAGRAM: 


SUBROUTINE RNPT 


Figure 38 
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and PRS2) and stress intensities (S1G1 and SIG2) in terms 
memune intersection coordinates, 5 and P. The values EM 
and BEE, which determine the pre-bucklinsg equation, are 


s x 


mpeed 1NCLO HPT via a COlMOR statement, along with other 


i] § pie 


normal program data inputs, including SIGA and SIGB. Input 
arguments include S1G1, SIG2 and PRS1, which along with the 
inelastic buckling pressure at point 2, PRS2 (computed by. 
RNPY from inout FCRE1), describe end points of the buckling 
2avectvion line approximation. 

Function BNLDS, This subprogram is the same as RNLDS 

Peqein the ideally plastic cass (i1.e., computation of elastic 

Lobar buckling pressure, PCRE1), with the addition of a 
promrammed iteration (figure 39, as used in reference 22) for 
Mme rnclastic portion. Zhe first Lwo stresses used in the 
iteration are SIGY and 0.9%SIGY. When the difference of 
computed inelastic failure pressure and assumed stress at 
which failure will occur in the next iteration becomes less 
em one half of one per cent of the present predicted 
feelure pressure, convergence is assumed. The eee ton 
sequence may be followed by using the diagram of the inelastic 
portion of RNLDS (figure 40) with the plot (figure 39). 

As noted before, identical procedures would be followed 
in the axisymmetric and peneral instability cases. The main 


: Ve aé yaa Serena : an lye is = To os =) a 4 > «. TINT ( 
Peer om would be unelCerea wit) tne exception, ayain ae notcd 
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ITERATIVE METHOD USED TO CONVERGE x42 


Clee otEC COLLAPSE PRESSURE... 









62,P) 


C 
ag) 





BUCKLING 
RQUATION 





6 
gee ty PREBUCKLING 
EQUILIBRIUM 
EQUATION 
cosy P.) 


COLLAPSE 
PRESSURE 


DIFFERENCE NEASURED 
APITER EAC Sena ttOnN 


Siz 09, (6, ,P,) 


COMPUTED INELASTIC FAILURE PRESSURES 
ASSUMED STRESS AT WHICH FAILURE WILL 
OCCUR IN THE NEXT ITERATION 


NOTE: METHOD TAKEN FROM REFERENCE 22 


Pees 99 
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(CONTINUING ON FROM 
PRESENT RNLDS PROGRAM) 1/42 


COMMON /G/ SIGA,SIGB /H/ BEE, EMH 


S1=0.9%S1GY 


CALI 
ROHOS 
Col eraM, boECs) 


[Pi=PCRP(PCRE1 ,T,ESEC,ETAN, SL) 


S2=SI1IGY 


CALL 
RNPT 
(BChinvensi, aol, 52, 
P2,P3,83) 


CAA 
Rie 
We Crim ol, 24.0 + 
S3,P4,P5,54) 














PCRP1=P4 





P3-PA> 0. 005uP3 NO. 


YES 
CALL 
RNPT 
(ecreing SL, Pes. 
Sue. P6,P7,S5) 


a 


P5-P6 > 0. 005#P 59 NO 


VES 
CALL 
RNPT 
Gael TE Sil,Po,o4, 
S5,P8,P9,86) 












P7-P8B > 0.005%P8 A NO 







neo S: 
CALL 
JGted ky 
GeC ie VS, roo), 
SOerLO. P11, S7} 


P9-P1050.005%P9% 







MGINLEY | 
4/10 Figure 40 





(CONTINUED FROM 










LAST PAGE...) me. 
'NON-CONVERGENCE! 
RNLDS=REDPR(PCRE1 ,PCRP1 ) 
(Figure 40 Continued) 
MCGINLEY 
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1h. 
Memeore, Of inse€rcion of Subprogram LINE. It is obvious 
that an optimization using the inelastic analysis would take 
much longer (verhavs by five or six times) than the ideally 


plasticcase. 
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